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ABSTRACT 

We present the analysis of phase-resolved X-ray and optical observations of the pecu- 
liar hot star HD 191612 (Of?p). This star is known to display line-profile variations 
that are recurrent with a period of 538 d and its spectrum was found to present the 
signature of a magnetic field. In the X-rays, it is slightly overluminous compared to the 
canonical Lx/ Lbol) relation and appears brighter when the optical lines are strongest. 
Our X.MM-Newton observations further reveal that the X-ray spectrum of HD 191612 
exhibits rather broad lines and is dominated by a 'cool' (0.2-0.6 keV) thermal com- 
ponent, two characteristics at odds with the proposed magnetic rotator model. We 
also report for the first time the low-level variability of the metallic (abs/em) lines 
and Hen absorptions that appear to be associated with radial- velocity shifts. Finally, 
we compare our results with observations of the early-type stars and discuss several 
possible scenarios. 

Key v^rords: stars: individual: HD 191612 - stars: early-type - X-rays: stars - X-rays: 
individual: HD 191612 



1 INTRODUCTION 

In 1972, a new subtype of massive stars, the Of?p category, 
was defined (Walborn 1972). It gathers stars that present 
C III emission lines around 4650 A of an intensity comparable 
to that of the neighbouring N ill lines. At that time, only two 
stars belonged to this category, HD 108 and HD 148937, but 
a third one was soon added, HD 191612 (Walborn 1973). 

Of the Of ?p stars, HD 108 has been the most studied 
in the past, but conflicting results were reported by differ- 
ent authors. To attempt to resolve the confusing situation, a 
long-term observing campaign was undertaken at the Haute- 
Provence Observatory (Naze et al. 2001). A detailed analysis 
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of the data allowed most of the older models to be discarded 
and showed that the star displays line-profile variations in 
the Balmer and the He I lines. Naze et al. (2001) also discov- 
ered that these line-profile variations seem to be recurrent 
with a timescale of approximately 50-60 years (see also Naze 
et al. 2006). 

Soon after, Walborn et al. (2003) reinvestigated an- 
other Of?p star, HD 191612, and discovered a similar phe- 
nomenon. Broad-band photometry of HD 191612 had been 
shown to exhibit a clear modulation with a period of 538 d 
(Koen & Eyer 2002; Naze 2004) and comparison with the 
spectroscopic results indicated that the spectral variations 
and photometric changes are correlated: when the star is 
faint, its emission lines are weakest and the star presents an 
apparent 08 spectral type, whereas when it is brighter, the 
emission lines are at maximum and the spectral type is 06.5 
(Walborn et al. 2004). 

Recently, Donati et al. (2006a) reported the detection 
of a magnetic field through the analysis of the Zeeman sig- 
nature in several lines of HD 191612, notably the variable 
He I lines. They proposed the star to be an evolved version 
of 6^ Ori C, the only O-type star for which a magnetic field 
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had been convincingly measured until then. Donati et al. 
(2006a) suggested that the 538 d period of HD 191612 was 
related to the rotation period of the star, which would have 
been braked as a result of the intense magnetic field. Un- 
fortunately, neither Donati et al. (2006a) nor other authors 
have yet observed a magnetic signature in HD 108 or mod- 
elled such a braking for this potential "twin" of HD 191612. 

To understand better the nature of Of?p stars, a de- 
tailed study of their high-energy properties is needed. For 
HD 108, observations were taken with XMM-A''eu)ton in Au- 
gust 2002, revealing a thermal X-ray spectrum (Naze et al. 
2004). The star also appears slightly overluminous in X-rays. 
Moreover, its X-ray luminosity was not significantly differ- 
ent from the previous detections by Einstein and ROSAT 
(noting, however, that the old data have very large error 
bars and that a decrease in flux by less than a factor of 2 
could not be excluded). These X-ray observations ruled out 
a compact-companion scenario where the wind ionization 
structure would be affected by the varying X-ray emission 
from an accretion disk (Naze et al. 2004). 

The shorter and better-known period of HD 191612 al- 
lows us to study its X-ray properties at crucial phases of the 
538 d cycle. Two years ago, a large observing campaign was 
therefore initiated to get multiwavelength, phase-resolved 
observations of HD 191612. We obtained four XMM-Newton 
pointings (PI: Naze) and many visible echelle spectra with 
the 1.93 m telescope of the Haute-Provence Observatory 
(PI: Rauw). We present in this paper an analysis of these 
two datasets. Additional optical data will be reported 
by Howarth et al. (in prep.), and throughout this pa- 
per, we wiU use their ephemeris: ro(HJD)=2 453 415.2±0.6, 
P=537.7±0.5d, where 0=0 corresponds to the maximum of 
the Ha emission (i.e. a phase shift of 0.5 compared to the 
ephemeris defined in Walborn et al. 2004). 

This paper is organized as follows. Section 2 describes 
the observations while Sections 3 and 4 present the results 
of the X-ray and optical OHP campaign, respectively. These 
observations are further discussed in Sect. 5 and we finally 
conclude in Sect. 6. 



2 OBSERVATIONS 

2.1 Haute-Provence Observatory 

Spectroscopic observations of HD 191612 in the visible 
domain were gathered during two years at the Haute- 
Provence Observatory (OHP). Most of the data were 
obtained with the 1.93 m telescope equipped with the 
Elodie echelle spectrograph, which covers the 3850-6800 A 
wavelength range at a resolving power of 42 000. Exposure 
times between 30 and 60 minutes were used, and the 
signal-to-noise ratio was generally about 100 or more. 
The Aurelie spectrograph at the 1.52 m telescope enabled 
us to get a few additional spectra. Three different set- 
tings were used: one covering the range 6360-6750 A at 
a resolving power R of 11000, the second 5480-5910 A 
at i?=9000, and the last 4450-4900 A at 7i=7000. Typ- 
ical exposure times were 15 minutes. All the reductions 
were performed using the MIDAS software developed at ESO. 



-4 







2 



_ , 1 , , , , 1 


|975 


1 


l9BI 


1 


• 

1004 




• 




/ 




1 






h 

|1068 # 



-0.5 0.5 





Figure 1. The phases of the XMyL- Newton observations (be- 
tween 2005 April and October) compared to the EW curve of the 
Ha line. Negative and positive EWs correspond to emission and 
absorption Hnes, respectively. The open circles indicate Aurelie 
data, filled symbols refer to Elodie observations, and crosses in- 
dicate other optical data (see Howarth et al., in preparation). 



2.2 XMM-Newton 

HD 191612 was observed four times with XMM-A'^emton for 
approximately 10-20 ks in 2005: Revolution (Rev.) 975 on 
Apr. 5, Rev. 981 on Apr. 17, Rev. 1004 on June 2 and 
Rev. 1068 on Oct. 8. For these observations, the three Eu- 
ropean Photon Imaging Cameras (EPICs) were operated in 
the standard, full-frame mode and a medium filter was used 
to reject optical light. The second observation (Rev. 981) 
was granted because the first data were affected by a rather 
long background flare; unfortunately, a long episode of high 
radiation took place at the time of this replacement obser- 
vation, rendering the data very noisy. 

We used the Science Analysis System (SAS) software, 
version 6.5, to reduce the X-ray data. They were first pro- 
cessed through the pipeline chains (tasks emproc, epproc 
and RGSPROC), and then filtered in order to keep only the 
best data: as recommended by the SAS team, we kept for the 
EPIC MOS (Metal Oxide Semi-conductor) detectors single, 
double, triple and quadruple events (i.e. pattern between 
and 12) that pass through the #XMMEA_EM fiffer; for the 
EPIC pn detector, only single and double events (i.e. pattern 
between and 4) with fiag=0 were considered. 

To check for contamination by low-energy protons, we 
have further examined the light curve at high energies 
(Pulse Invariant channel numb er> 10000, E>10keV, and 
with pattern=0). Where background flares were detected, 
we discarded the time intervals presenting high count rates 
(larger than 0.6 cts s"^ for EPIC MOS, 4.0 cts s"^ for EPIC 
pn, and 4-7.5 cts s~^ over the full energy range of the Re- 
flection Grating Spectrometer, RGS) . The best observations 
(i.e. Revs. 975, 1004, and 1068) were also combined with 
the tasks merge for EPIC or rgscombine for RGS. Fur- 
ther analysis was performed using the SAS, the FTOOLS 
tasks and the XSPEC software v 11.2.0. A table, available 
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electronically at the CDS , provides the X-ray properties 
of all the sources detected in the field. 



3 X-RAY RESULTS 

3.1 Light-curve 

The X-ray observations were planned to sample crucial 
phases in the 538 d period (Fig. 1): the first observation was 
taken near the time of maximum emission in the visible do- 
main while the last one was made during the 'quiescent', 
low emission phase. All the X-ray data were taken during 
the same 538 d cycle. 

The count rate of HD 191612, evaluated using the task 
EDETECT_CHAIN, decreased by 40% between the first (Apr. 
5, Rev. 975, (f>^0.1) and the last (Oct. 8, Rev. 1068, cf>^0.4) 
observations, as shown In Fig. 2. At the same time, simulta- 
neous visible observations taken at the OHP display a weak- 
ening of the emissions in the H and He i lines, as expected. 
Prom a detailed study of the behaviour in the visible of 
HD 191612 (Walborn et al. 2004), wc know that these spec- 
troscopic changes correspond to a decrease of the broad- 
band luminosity. Within the limitations of our dataset, our 
observations therefore suggest that the X-rays follow the 
behaviour observed in the visible range. Further data are 
however needed to fully assert this conclusion. 

Note that such a modest decrease in the X-rays could 
not be detected in the low signal-to-noise ROSAT data 
of Of?p stars; the result is thus consistent with the non- 
detection of such variations in HD 108 (Naze et al. 2004). 
Finally, on short timescales (i.e. during one pointing), we de- 
tected no significant variations of the X-ray emission. In fact, 
O stars generally lack variations on timescales this short 
(see, e.g., Sana et al. 2006). 

Looking at the hardness ratios of HD 191612 (Fig. 2), 
two different behaviours appear, depending on the energy 
range considered: at lower energies (0.4-2.0 keV), the de- 
creasing X-ray emission clearly softens with time whereas 
the proportion of the harder emissions in the 1-2 keV and 
2-10 keV ranges does not change significantly with phase. 

3.2 Low-Dispersion Spectroscopy 

EPIC spectra of HD 191612 were extracted over a circular 
region with radius 30" centred on the star. A faint source 
is present near HD 191612, so we did not choose an annu- 
lar background region, but rather a nearby circle devoid of 
sources. For a colour excess E{B — V) of 0.59 (Savage et al. 
1985), the interstellar absorbing column can be estimated 
as A^i^t = 3.4 X 10^^ cm~^ (using the relations of Bohlin 
et al. 1978). In the spectral modelling, we do not allow the 
absorption to go below this threshold. 

The spectra from each instrument were fitted sepa- 
rately, but as they gave consistent results, within the er- 
rors, we finally fit all EPIC data simultaneously (Fig. 3). 
Fits using only one thermal component were unacceptable 
(x^ > 2). We rather used the sum of two absorbed, optically 
thin equilibrium plasma models {mekal, Kaastra 1992). 
Since these two thermal components could arise in different 




Figure 2. Variations with phase of the count rate in the 0.4- 
10.0 keV band (top) and of two hardness ratios (middle, bottom). 
Open symbols (triangles and squares) refer to EPIC MOS data, 
whereas filled circles represent EPIC pn observations. The hard- 
ness ratio HRi (resp. HR2) is defined as {M-S)/{M + S) (resp. 
{H - M)/{H + M)), where S is the count rate in the 0.4-1.0 keV 
band, M in 1-2 keV, and H in 2-10 keV. 



regions of the wind, we allow them to be absorbed by inde- 
pendent column densities (N^). The best fit was obtained 
when the cooler component dominates, as is typically ob- 
served for 'normal' 0-type stars (see, e.g., Sana et al. 2006). 
In fact, two sets of parameters gave similar values of the 
chi-squared (see Table 1), one with temperatures kT around 
0.6 and 1.0-2.5 keV and negligible absorbing columns N^; 
the other with lower temperatures (0.24 and 1.0-1.6 keV) 
and larger columns {N^ ~ 0.5 and 1 x 10^^ cm"^). The lat- 
ter solution presents slightly lower x^i even when including 
RGS data. Moreover, a fit of the EPIC data with a differen- 
tial emission-measure model (DEM, cdpmekl, Lemen et al. 
1989; Singh et al. 1996) confirms the existence of a domi- 
nant temperature near 0.2 keV whilst a secondary hump in 
the DEM distribution appears at 1-2 keV (Fig. 3). 

Using bolometric corrections from Humphreys & McEl- 
roy (1984) and the optical properties of HD 108 {V=7A, 
B - V=0.18, d=2.51 kpc) and of HD 191612 (1^=7.8, B - 
V=0.27, d=2.29 kpc, Gies 1987), we calculated the bolo- 
metric luminosities of both objects: I/bol=1.6x 10"^^ erg s^^ 
for HD108 and 1.2x10^'' erg s"^ for HD 191612. The X- 
ray luminosities Lx"°''"*, evaluated from the EPIC data in 
the 0.5 -10.0 keV range and corrected for interstellar absorp- 
tion, are 1.2x10^^ erg s"^ for HD 108 (Naze et al. 2004) 
and 7 - 9 X 10^^ erg s"^ for HD 191612 (this work). This 
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results in log(L^"*7LBOL) ratios of -6.1 for HD 108 and 
-6.1 to -6.2 for HD 191612. The values of the L'^'^^yLBOh 
ratio arc thus five times larger than that of the 'canonical' 
relation (log(Z/x"''''V-^'BOL) = -6.91 in the most recent de- 
termination of Sana et al. 2006), as shown in Fig. 4. This 
over luminosity is slight but significant, since the new study 
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HD191612: EPIC data from Rev1068 HD1S31612: KPTC data from Rev. 1068 
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Figure 3. Le/t.X-ray spectra of HD 191612 recorded by EPIC-pn (upper data) and MOS (lower data) during Rev. 1068 (</)=0.4, Oct. 
2005), with the best- fit model superimposed (see Table 1 for values of the parameters). Right: Differential emission measure distribution 
from a fit of the same EPIC data by a DEM model where <j>{T) corresponds to the exponential of a sum of Chebyshev polynomials 
between the 1st and 6th orders. The error bars were derived from Monte-Carlo simulations. 
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Figure 4. Diagram showing the X-ray luminosity (in ergs s~^) 
versus bolometric luminosity (in ergs s~^). The dashed line in- 
dicates the typical relation for O stars (from Sana et al. 2006), 
which HD 108, HD 191612 and 6»l Ori C all lie above. Asterisks 
show the position of hot stars in NGC 6231 (Sana et al. 2006) with 
three outliers: the two objects lying above the line are coUiding- 
wind binaries whereas the one lying below is a Wolf-Rayet binary. 



by Sana et al. (2006) shows that the scatter around the 
'canonical' relation is only about 20%. 

The parameters found for HD 191612 (two tempera- 
tures, soft/hard character, overluminosity) are very similar 
to those found for HD 108 by Naze et al. (2004), confirming 
the similarities between these two Of?p stars, even in the 
X-ray domain. 
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Figure 5. Smoothed combined RGS spectrum of HD 191612 be- 
tween 8 and 20 A (RGSl -I- RGS2, first and second order. Revs. 
975, 1004, and 1068). The strongest lines are identified. 



3.3 High-Resolution Spectroscopy 

Due to the weakness of HD 191612 in X-rays, individual RGS 
spectra do not have enough counts for a sound analysis (see, 
however, below). We have therefore first combined the data 
from the three observations with the lowest background lev- 
els (i.e.. Revs. 975, 1004, and 1068) using the task rgscom- 
BINE and then assembled the first and second-order spectra 
of both instruments into a single RGS spectrum by means 
of the RGSFLUXER command. 

Figure 5 shows the resulting RGS spectrum in the wave- 
length range 8-20 A. Due to the rather heavy interstellar and 
circumstellar absorption, and the weakness of the source, the 
signal-to- noise ratio of RGS data above ~ 20 A is too low 
to be useful. The spectrum of HD 191612 is clearly domi- 
nated by emission lines and, apart from some higher excita- 
tion lines (e.g. Fexvill) that appear relatively strong here, 
it looks quite similar to the RGS spectrum of typical O stars 
such as C, Pup (Kahn et al. 2001). In particular, we note that 
the RGS spectra do not show a strong bremsstrahlung con- 
tinuum emission unlike what is seen in the HETG spectra 
of the colliding-wind binary WR140 (Pollock et al. 2005) or 
the peculiar mag netic O star Ori C (Schulz et al. 2000). 

We identified the strongest lines by comparison with the 
line list of the SPEX plasma code (Kaastra et al. 2002) . The 
most prominent lines are hydrogen- and helium-like transi- 
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Table 1. Best-fitting models and X-ray fluxes at Earth for each XMM-Afe«;4on observation of HD 191612. The fitted model has the form 
wabs(AfPj.)*[wabs(Af{l)*mekal(kTi)+wabs(7V|l)*mekal(kT2)] , with wabs(iVPj.)= 3.4 X 10^^ cm"^. Quoted fluxes are in the 0.4-lO.OkeV 
energy range. The unabsorbed fluxes f^^^^ are corrected only for the interstellar absorbing column. For each parameter, the lower 
and upper limits of the 90% confidence interval (derived from the ERROR command under XSPEC) axe noted as indices and exponents, 
respectively. The normalisation factors are defined as ^'^^^ / neUndV, where D, Ue and nn are respectively the distance to the source, 
the electron and proton density of the emitting plasma; (p is the phase in the 538 d cycle. 



Date 


Rev. 




10^2 cm-2 


kTi 
keV 


normi 
lO-^cm-S 


10^2 cm-2 


keV 


10-3cm-5 


xS(dof) 


(10-13 


j?unabs 
ergcm-2 s-l) 


'cool' model 






















05/04/05 
17/04/05 

02/06/05 
08/10/05 


975 
981 

1004 
1068 


0.09 
0.12 

0.20 
0.44 


U.OIq 42 

0.32«:« 

o.43[;:i; 

e,^o■6'i 

^•"^0.58 


U.ZOq 20 

n 27" -29 

n 250-28 
'J-^'JO.23 

o.2o8:?i 


7 0715-9 
'•'J'4.04 

2 1 81 
3-'Ol.78 

i3.6ig:i 


1 07I-31 
-'-•'J'0.87 

qil-io 

1 031-25 

1 QQl.68 

i.jyj ^4 


1 27I-3'' 
-'-•^'1.19 

1 0ll-"7 

1 991.31 

1 60l-'i'2 
LDU]^ 49 


-1 991.35 
-■-•^^1.11 

0.891:11 

0.668:5! 


1.16 (411) 
1.10 (681) 
1.15 (460) 
1.31 (555) 


7.0 
6.8 
6.0 
5.3 


14.1 
14.3 
12.5 
11.0 


'hot' model 






















05/04/05 
17/04/05 
02/06/05 
08/10/05 


975 

981 
1004 
1068 


0.09 
0.12 
0.20 
0.44 


n 0.05 

n 0.01 
'J-0. 
n 0.01 
^•0. 

0.01 


0.648:66 

n fin^-oi 

"■'^'^0.58 
n cil0-54 

o-6i8:ig 


0.288:33 

o.328:if 
o.288:i 

n 9'<0.24 
"■^•^0.22 


n 0.03 

0.818:70 

0.828:fi 

0.8:°! 


9 Qc;2.B4 

-1 90I.32 
^^■^■^1.06 
1 07^-27 

J^-^'i.oi 

9 c;o2.68 


o.5i8:ti 

o-95i:ii 

921-01 
^•^^0.75 

n Q9O.34 


1.19 (411) 
1.08 (681) 
1.16 (460) 
1.39 (555) 


7.6 

7.0 
5.8 
5.4 


14.4 
13.9 
12.0 
11.0 



NelX FeXVIII FeXVII 




Revl068 
= 0.4 



X (A) 



Figure 6. Smoothed combined RGS spectra of HD 191612 (RGSl 
+ RGS2, first and second order). The spectra from Rev. 1004 
and 975 are shifted upwards by 7 x 10-^ and 1.4 x 10-* units 
respectively. Note the changing intensity of the Fexvill lines at 
AA 11.421, 14.204, 14.212 and Nex Lya compared to the neigh- 
bouring lines of Neix (helium-like triplet) and Fexvii. 



tions from Mg, Ne, and O together with hnes from FeXVII 
and Fe XVIII (see Fig. 5). We also note the possible presence 
of weak spectral signatures of additional iron ions; for ex- 
ample, the strongest Fexix line (A 13.52) might be present, 
but would be heavily blended with the strong Ne ix He-Uke 
triplet. The strongest Fexx lines could actually be responsi- 
ble for the low-level features near 12.8 A, whilst the strongest 
Fexxi line (A 12.29) could account for some of the breadth 
of the Nex Lya line. The spectral signature of these ad- 
ditional ionization stages of iron is indeed expected since 
the EPIC spectra reveal the likely presence of an Fe XXV 
lino at 6.7 koV. The temperatures of maximum omissivity of 
the most prominent lines seen in the RGS spectra span a 
range from 2 - 3 x 10^ K (O vii and O viii) to 7 - 10 x I06 K 
(FeXVIII and Mgxil). These values, in agreement with the 
results from the low-dispersion spectroscopy, arc again very 
similar to those found for the Of?p star HD 108 (Naze et al. 
2004). 

The rather limited signal-to-noise ratio of the spectrum 
prevents us from doing a detailed profile analysis for in- 



dividual lines. We have nevertheless measured the fwhm 
of Ovill Lya and the strong Fexvil A 15.014 line on the 
combined (unsmoothed) RGS spectrum, finding fwhm val- 
ues between 2000 and 2600kms-i. We also tentatively per- 
formed a global fit of the X-ray lines by a single triangular'^ 
profile whose amplitude was allowed to vary from line to 
line but whose width and shift were taken to be the same 
for all lines. Within the limit of our RGS data, we detect no 
blue/red shift >250kms^inor any significant line asymme- 
try. The FWHM is found to be about 1800±400 kms'i, or 
70-85% of the terminal velocity (Walborn et al. 2003), sug- 
gesting again rather broad lines. As the above values corre- 
spond to about twice the instrumental resolution, it seems 
that the spectral lines of HD 191612 are indeed broad and 
can be resolved with the RGS. 

The strength of some of the spectral lines changes be- 
tween the various XMM-Newton pointings. For instance, the 
Fexviii lines appear strongest in the spectrum from June 
2 (Rev. 1004, 0=0.2), when the star is in an intermediate 
state; at that time, Fexviii AA 14.212, 14.260 have an inten- 
sity comparable to that of Fexvii A 15.014. However, these 
lines are weaker in the first and last observations (Revs. 975 
and 1068, 0=0.1 and 0.4; Fig. 6). A similar difference is seen 
for the Ne x and Ne ix lines. These variations most probably 
reflect the changing hardness ratio measured from the EPIC 
data of HD 191612. 



4 PROPERTIES IN THE OPTICAL RANGE 

To identify the lines displaying variability, we computed a 
Temporal Variance Spectrum (TVS, Fullerton et al. 1996) 
for the Elodie data set (shown partly on Fig. 7) . As already 
mentioned in Naze et al. (2001) for HD 108 and in Walborn 



2 The triangular profile gives an unbiased estimate of any shift 
and asymmetry and is similar to several proposed theoretical pro- 
files. 
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Figure 7. Mean spectrum of HD 191612 (top panels) and the TVS (bottom panels). The main lines are indicated at their rest wavelengths. 

For the TVS panels, a dotted horizontal line indicates the significance level at 99%. Some variability due to changing weather conditions 
appears at the position of the telluric lines, e.g. near Cm A 5696. 
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Figure 8. Variations of HQ (left panel) and C IV AA 5801,5812 
lines (right panel). The individual panels display, from top to 
bottom, the mean line profile, the TVS, as well as the amplitude 
and the phase of the sine wave of frequency 1/537.7 d'^ fitted to 
our data. 



et al. (2003) for HD 191612, the most significant profile vari- 
ations are found in the Balmer and Hei lines. The TVS of 
these lines displays a roughly gaussian shape centered on 
the rest wavelength. It is important to note that the blue 
wings of the profiles, unaffected by the emission, exhibit 
very little variability (Fig. 8). The Balmcr and He I lines 
thus consist of an appaxently constant^ absorption compo- 
nent with a superimposed varying emission profile extending 
from approximately —200 to 200kms~^. For 0-type stars, a 
strong C III AA 4649,4650 emission is quite rare. Its presence 
in HD 191612 cannot be due to an abundance effect (Wal- 
born 2001) because the strength of these lines varies while 



^ Because of the presence of the varying emission, small RV vari- 
ations as those exhibited by the metallic lines (see below) would 
not be detectable for the observed H and He I lines. 



C in A 5696 remains constant. Since the emission mechanism 
of C III AA 4649,4650 is rather complex (Andrillat & Vreux 
1975), the explanation of their variability requires detailed 
atmospheric modelling. 

On the other hand, now that we have data with high 
signal-to-noise covering the whole 538 d period, we see that 
almost all lines present significant changes (Fig. 7) . The vari- 
ability characteristics, though, are different from those of the 
Balmer and He I lines since their TVS generally displays a 
double-peaked shape (Fig. 8). 

We have analyzed further the spectral variations by us- 
ing techniques similar to those of Rauw et al. (2001). At ev- 
ery wavelength stop, we have fitted a sine wave of frequency 
1/537.7 d~^ to our time series. Fig. 8 shows the amplitude 
and the phase of this sine wave as a function of radial ve- 
locity for two different spectral regions. Note that the phase 
(j)o has no physical meaning outside the wavelength domain 
where a significant variation is detected. This explains the 
noisy phases appearing at the extremes of our velocity in- 
terval. 

Once again, the two subsets of lines display different be- 
haviours, in phase as well as in amplitude. Absorption lines 
from Hen and abs./em. lines from metals present a double- 
peaked amplitude profile and a varying phase, whereas for 
the Balmer and Hei lines, the amplitude profile is gaussian 
and the phase is rather constant (Fig. 8). We note however 
that a small double-peaked feature is apparently superim- 
posed on the gaussian amplitude profile for the weakest He I 
lines and He ii A 4686. 

To further study the spectral behaviour of HD 191612, 
we decided to focus on the C iv AA 5801,5812 lines since (1) in 
our data, these lines belong to the domain presenting the 
highest signal-to-noise ratio and (2) they display a very clean 
gaussian shape that does not seem affected by any resid- 
ual emission or distortion of the line profile. These lines arc 
thus most likely to arise in the stellar photosphere, unlike 
the Balmer and Hei lines which are often contaminated by 
variable wind emission. The results obtained from these C IV 
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Figure 9. Left panel: Radial velocity (top) and equivalent width 
(bottom) measurements for tlie C IV A 5801 line. Filled and open 
symbols stand for Elodie and Aurelie data, respectively. Right 

panel: Same for He II A 4686. The RV points corresponds to the 
highest emission peak. Typical errors are 1— 2kms^^ for RVs and 
0.05A for EWs. 



lines should therefore remain unaffected by strong noise or 
changing wind properties. 

First, the equivalent widths (EWs) were measured by 
integrating the line profiles and we found that they remain 
rather stable, with an apparently stochastic variability of 
maximum 10%. A similar result was found for other prob- 
able photospheric lines (i.e. metallic lines and He ii absorp- 
tions). In marked contrast, the variable hydrogen. He I and 
He II A 4686 lines follow a very different behaviour, as was de- 
scribed in Walborn et al. (2004) for Ha: their EWs first rise 
smoothly, reach a maximum emission level and then decline 
smoothly (Fig. 1). 

The radial velocities (RVs) wore determined by fitting a 
gaussian to the C iv line profile. They appear systematically 
larger in the last OHP observations (Fig. 9). Again, similar 
results were found for the other strong metallic lines and 
the He ii AA 4542,5412 absorption lines (where the RVs were 
derived from a gaussian fit to the bottom - resp. top - of the 
absorption - resp. emission - profile). 

Since the C iv lines are most probably of photospheric 
nature, their apparent motion cam a priori be attributed to 
the star itself with some confidence. Although the periodic 
behaviour of these RVs cannot be assessed here with confi- 
dence, such a trend in RVs together with the double-peaked 
profile of their TVS suggest a binary nature for HD 191612 
(although probably with a period differing from 538 d). 

The case of the Hell A 4686 emission line is peculiar: 
the RV of its emission peak follows the same trend as the 
RVs of metals when the emission is minimum (Fig. 9). This 
suggests that there are two components to the emission of 
He II A 4686, one varies in harmony with the Balmer lines 
(and thus represents an additional omission) while the other 
simply shifts in wavelength as for the metallic lines (thus 
probably belonging to the star). Note that an additional 
component, in absorption, is also clearly detected. This 
might indicate the presence of a companion. 



5 DISCUSSION 

After presenting the detailed analysis of the pliasc-rosolvod 
observations taken in the X-ray and optical ranges, we now 
discuss the implications of these results on the possible na- 
ture of HD 191612. 



5.1 HD 191612, a magnetic oblique rotator? 

Because of the detection of a magnetic field, one may think 
of HD 191612 as an analogue to the other two magnetic hot 
stars 6^ Ori C and r Sco (see e.g. Donati et al. 2006a). In- 
deed, these systems share several similarities (see Table 2 
for a comparison). For example, they display large magnetic 

confinement parameters {rj = -^j-^) and their Alfvcn radii 
are also similar (about 2 R, ; Donati et al. 2006a,b). In addi- 
tion, their log(iJ"^''7LBOL) ratios, -6.0 for 9^ Ori C and up 
to —6.1 for HD 191612, indicate an important overluminos- 
ity in X-rays compared to 'typical' O stars (Fig. 4). Finally, 
these three magnetic objects display variations in the X-rays 
that seem correlated to those observed in the visible domain 
(Gagnc et al. 2005, this work). 

However, our new data also reveal many differences. 
First of all, the X-ray spectra of 9^ Ori C and r Sco are 
much harder than that of HD 191612. The temperatures 
of the hot plasma of 9^ Ori C amount to 0.7 and 2.6 keV 
(i.e. 8 and 30MK; Gagne et al. 2005) where the higher- 
temperature component clearly dominates. The opposite sit- 
uation is found for HD 191612, even in the case of the 'hot' 
fit. This distinct behaviour cannot be attributed to extrinsic 
causes since both systems share similar interstellar absorp- 
tion: the hard X-ray emission of HD 191612 is thus unlikely 
to be hidden by a high extrinsic column density. 

Differences are also evident when considering r Sco, 
which has a somewhat softer spectrum than 0^ Ori C. In 
this case, the differences can be best appreciated by com- 
paring the differential emission-measure fit given by Mewe 
et al. (2003) with ours: again, the hard component is much 
reduced for HD 191612. The 'hard' character of the X-ray 
emission of 9^ Ori C and t Sco is also revealed by the pres- 
ence of several strong, high-ionization lines (e.g. Fexxiv) 
which are weak or absent in the spectrum of HD 191612. 

We also note that 9^ Ori C presents a prominent 
bremsstrahlung continuum in its X-ray spectrum which is 
a clear indication that the emission-measure distribution is 
dominated by high temperatures. No such feature is seen in 
the case of HD 191612. 

In addition, the X-ray emission of 9^ Ori C and 
T Sco consists of very narrow lines (FWHM well below 
1000 km s"^ see Gagne et al. 2005; Cohen et al. 2003), as 
predicted theoretically for a magnetically confined wind, but 
this is not the case for HD 191612, as was shown in Sect. 3.3. 
Even the broadest lines of 9^ Ori C and r Sco (notably O vill 
Lya and the strong Fexvil A 15.014 line) are still narrower 
than the lines of HD 191612. The line width measured for 
HD 191612 is actually much more similar to that determined 
for 'normal' O stars. In many respects, HD 191612 thus ap- 
pears as an intermediate case between the magnetic objects 
and the typical O stars (Fig. 10). 

Finally, we note that the behaviour of HD 191612 in the 
visible range also differs somewhat from that of 9^ Ori C, al- 
though it is difficult to compare them in detail. In fact, while 
the spectral variability of HD 191612 is now well known, 
this is not always the case for 9^ Ori C: for example, vari- 
ations of the He I A 4471 /He ii A 4542 ratio of 9^ Ori C were 
reported by Walborn (1981) but the same ratio was subse- 
quently found constant (Stahl et al. 1993). 

The differences between HD 191612 and Ori C or t 
Sco, especially in the X-ray domain, are quite puzzling, in 
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Table 2. Properties of 9^ Ori C (from Gagnc ct al. 2005), HD 191612 (from Walborn ct al. 2003; Donati et al. 2006a), and 
rSco (from Donati et al. 2006b). The wind density at i'tx>/2 was calculated assuming a /3-type law (with ^=0.8) and 
is shown in units of 10"'^'* g cm~^. The X-ray luminosity corresponds to the 0.5-10.0 keV band. 
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R, 




M 






Pw 




V 


r unabs 


Period 






(Re) 


(kK) 


(Mq yr - 




(km s~l) 




(G) 




(erg s-i) 


(d) 


HD 191612 


06 


16.7 


37 


6.1x10" 


-6 


2600 


5.9 


-1500 


7.6 


7-9x10^2 


538 




08 


18.7 


35 


2.5x10" 


-6 


2100 


2.4 




29. 






el Ori c 


05.5V hot 


8.3 


45 


1.4x10" 


-6 


2980 


4.7 


1060 


3.6 


0.8-1.1x10^^ 


15 




05.5V cool 


9.1 


42 


5.5x10" 


-7 


2760 


1.7 




12 






T Sco 


BOV 


5.2 


32 


2x10" 


8 


2000 


0.3 


not dip. 


~40 


2-3x10^1 


41 



view of other similarities. One could imagine that the larger 
mass-loss rate of HD 191612 compared to 9^ Ori C would 
prevent us from seeing the X-ray emission from material in 
the magnetosphere (i.e. below the Alfven radius). However, 
recent analyses of the X-ray line profiles of (" Pup (Kran- 
mer ct al. 2003) and C, Ori (Cohen et al. 2006) indicate that 
in these stars (which are more similar to HD 191612 as far 
as M is concerned). X-ray emission from down to 1.5 R» 
can be observed. Therefore, if there is a magnetically con- 
fined wind region in HD 191612, we would expect to be able 
to see its X-ray emission, unless the actual magnetic field 
of HD 191612 is much lower than the dipolar one (1.5kG) 
inferred by Donati et al. (2006a). This is indeed possible 
since the exact magnetic configuration of HD 191612 is not 
yet known (only one phase was observed by Donati ct al.). 
The discrepancies found between the actual properties of 
HD 191612 and the magnetic oblique rotator model (as rep- 
resented by 6^ Ori C) indicate that the emission mechanisms 
of these peculiar stars is not fully understood yet. 



inclination system with a lower-mass companion. Owing to 
the youth of HD 191612, this companion, if stellar, should 
still be pre-main-sequence (PMS) but the low luminosity of 
such objects {My ~ -2.5 for a TM© star, Siess et al. 2000) 
compared to mid-O-type stars could prevent its emission 
being seen in the composite spectrum of the system. The 
alternative case of an accreting, compact object seems more 
promising. Of course, high accretion rates can be readily 
ruled out since they would lead to a much larger X-ray lu- 
minosity (see Naze et al. 2004 and references therein) . How- 
ever, the observed slight overluminosity in X-rays might be 
compatible with low to moderate accretion rates. Variations 
of the optical and X-ray emissions would be explained by the 
eccentricity of the orbit. However, this model also encoun- 
ters difficulties when analyzed in detail: for example, such 
an eccentric orbit would result in asymmetric EW curves for 
Ha and similar lines, which is not the case. 



5.2 Alternative Scenarios 

Since there are still problems with the magnetic oblique ro- 
tator scenario, we might envisage alternative models. If we 
assume that the RV variations are periodic and indeed pro- 
vide the signature of a binary, the low amplitude of the RV 
curve (30kms~^) would suggest a low mgiss function, which 
leaves us with two possibilities. 

First, that the system of HD 191612 is seen at a low 
inclination and is then composed of two 'normal' early-type 
stars. In this case, one may immediately think of coUiding- 
wind (CW) phenomena. Indeed, Sana et al. (2006) found 
that large log(Lx"'''"*/-^BOL) ratios were displayed by the 
only two CW systems in their sample (sec Fig. 4). In ad- 
dition, the presence of a high-temperature component in 
the X-ray spectrum, and the modulations of the optical and 
X-ray emissions, are generally considered as typical signa/- 
tures of this phenomenon. However, a detailed hydrodynam- 
ical simulation would be required to predict the putative 
coUiding-wind emissions and their variations; while further 
X-ray observations are necessary to define better the cor- 
relations between X-ray and optical variations. In addition, 
no obvious signature of a massive companion has yet been 
detected in optical observations, except maybe for the pe- 
culiar shape of the He ii A 4686 profile (but see Howarth et 
al., in prep.). 

On the other hand, HD 191612 could be a high- 



6 CONCLUSIONS 

We have obtained phase-resolved observations of HD 191612 
with XMM.-Newton and the Elodie/Aurelie spectrographs 
of the OHP. 

In the high-energy domain, the star presents the spec- 
trum of thermal plasmas with temperatures of 0.2-0.6 and 
1-2 keV. The cooler component dominates, and no strong 
lines with ionization higher than FeXVIII are observed. The 
width of the X-ray lines, measured in the high-resolution 
RGS spectra, is 1800±400 km s"^ These results, which are 
similar to those found for 'normal' O stars, are at odds 
with the predictions of a simple magnetic rotator model. 
Although no variability is detected during each pointing, 
the X-ray flux decreased by 40% between the first and last 
XMM-iVewton observations which were taken during high 
and low emission states (as defined in the visible), respec- 
tively. The high-energy variations might thus follow the op- 
tical changes. The average X-ray flux exceeds by a factor of 
5 that derived from the 'canonical' L'^'^^^ / Lboi. relation. 

A clear dichotomy is seen in the behaviour of the optical 
lines. On the one hand, the hydrogen and He I lines display a 
relatively constant absorption superimposed on an emission 
of varying intensity (centred on the rest wavelength). On the 
other hand, our high quality dataset enabled us to detect 
for the first time the variations of the He ii absorptions and 
metallic (abs. or em.) lines. This variability is apparently 
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Figure 10. EPIC MOS spectra of C Pup, HD 191612, and 
0^ Ori C: HD 191612 appears as an intermediate case between the 
"typical" O-type star (soft spectrum) and the oblique magnetic 
rotator (hard spectrum). The spectra have been arbitrarily re- 
normalized to appear with similar strengths, and an extra absorp- 
tion was added to have all stars under AfH=4.5xlO^^ cm"'^, there- 
fore correcting for the dissimilar interstellar absorptions of these 
objects (log{Nu) of 20.0, 21.5, and 21.7 for C Pup, HD 191612, 
and 9^ Ori C, respectively). 

associated with shifts in radial velocity (up to 30kms^^) 
which are reminiscent of a binary motion. 

The He ll A 4686 line displays a mixed behaviour which 
could tentatively be explained by the presence of three com- 
ponents: an emission following the behaviour of the H and 
He I lines, an emission shifting in RV as the metallic lines, 
and one additional absorption, probably of stellar origin. 

With these new data, the mystery deepens on the nature 
of Of?p stars: HD 108 and HD 191612 display some proper- 
ties of typical O-type stars but also share some similari- 
ties with peculiar objects such as 6^ Ori C (Fig. 10). The 
puzzling results from our multiwavelength campaign thus 
challenge the models proposed to explain the behaviour of 
HD 191612. 
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